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Selected  Studies  of  Underwater  Sound  Phenomena 


Preface 


/  This- document  reports  on  three  separate  studies  conducted  during 
fiscal  year  1966.  Pagination,  references  and  figures  are  numbered 
separately  for  each  study. 


Part  I,  Some  .Stochastic  Structures  for  Modeling  Underwater  Sound, 
discusses  simulation  of  propagation,  reverberation  and  ambient  noise 
with  an  underlying  stochastic  description  of  the  sea  surface. 

Part  II,  Spectra  of /The  Dynamic  Sea  Surface  proposes  a  unified 
model  of  the  equilibrium  sea  surface,  including  turbulence  and  based  on 
observations  of  several  different  features  of  the  sea. 


Part  m,  System  Parameters  for  Ultrasonic  Explosive  Echo-Ranging, 
illustrates  the  significance  of  high  frequencies  (up  to  40  kHz)  to  permit 
improved  directional  discrimination. 


Some  Stochastic  Structures  for  Modeling 
Underwater  Sound 


H.  W.  Marsh* 


Abstract 

Seversil  stochastic  structures  are  described  which  may  be 
used  to  model  certain  features  of  the  underwater  sound  field. 
These  include  propagation  { and  reverberation)  and  sea  noise. 
Emphasis  is  placed  upon  models  which  are  adapted  to  digital 
computation,  in  which  wave  forms  are  of  interest,  in  addition 
to  averages  such  as  power  spectra. 
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Introduction 


The  principal  source  of  stochastic  effects  upon  underwater  sound  is 
the  dynamic  sea  surface.  There  are,  of  course,  other  water  motions,  in¬ 
cluding  tides,  internal  waves  and  currents.  However,  these  generally  tend 
to  produce  rather  smaller  effects  and  of  longer  periods  than  the  wind  driven 
sea.  Furthermore,  in  a  compromise  between  scope  and  detail,  our  atten¬ 
tion  is  going  to  be  focused  on  a  simple,  but  important,  underwater  sound 
situation  in  which  the  presence  of  the  sea  bottom  and  of  sound  velocity 
structure  are  not  emphasised.  We  shall  not  treat  sea  noise  other  than  that 
which  is  generated  by  the  sea  surface.  Under  these  conditions,  a  single 
stochastic  function,  the  turbulence  spectrum  of  wind  driven  waves,  is  suf¬ 
ficient  to  permit  modeling  of  noise,  propagation  and  reverberation.  We 
shall  discuss  the  turbulence  spectrum,  and  show  a  model  of  the  sound  field 
which  can  be  expected  to  produce  covariance  (space-time  )  functions  which 
are  indistinguishable  from  measurements.  More  than  this  is  obtainable 
from  the  model:  wave  forms.  While  the  wave  forms  cannot  be  defended 
upon  empirical  grounds,  since  measurements  are  very  limited,  they  will 
at  least  serve  as  a  guide  to  what  may  be  expected.  Further,  by  following 
in  the  footsteps  of  Stocklin1  ,  the  field  continuum  may  be  recovered  to  the 
extent  shown  by  his  theory  of  soace-time  sampling.  The  modeling  will  be 
described  in  a  hybrid  context  This  permits  interspersed  use  of  actual 
sound  recordings  and  digital  simulation,  as  well  as  flexibility  in  choosing 
specific  mechanisations.  We  have  not  yet  produced  an  entire  simulation, 
but  will  show  examples  of  results  to  date. 


The  Physical  Situation 


A  simple  found  source  is  located  at  the  origin  of  coordinates.  At  a 
distance  R  from  the  source,  the  free  field  pressure  f  (t  -  R/c  )/R  is 
assumed,  for  constant  sound  speed  c.  If  non-linear  propagation  effects 
can  be  neglected,  an  extended,  directional  source  may  be  considered,  and 
near  field  effects  included  without  difficulty.  However,  in  order  to  accom¬ 
modate  non-linear  effects,  we  shall  deal  here  with  source  characteristics 
dependent  only  upon  the  radial  coordinate. 

At  any  point  in  the  medium,  with  position  vector  R  (x,  y,  r)  the  total 
sound  field  will  be 

p(R,  t)  =  Pa  +P d  +P#  (1) 

P  _  denotes  the  ambient  pressure,  entailing  the  hydrostatic  head  as  modu¬ 
lated  by  the  moving  sea  surface.  and  P#  are  produced  by  the  source, 
and  denote  respectively  a  direct  arrival  and  an  arrival  scattered  by  the  sea 
surface.  P#  in  turn  will  consist  of  a  coherent  component  P  corre  - 
s ponding  to  ordinary  reflection  with  inversion  of  sign  and  loss  of  amplitude, 
and  an  incoherent  component  P  .  The  latter  is  "reverberation". 

The  sea  surface  lies  at  a  mean  elevation  Zq  above  the  sound  source 
and  is  described  by  the  instantaneous  elevation  Z  (x,  y,  t ).  Its  form  and 
motion  determines  the  ambient  pressure  and  scattered  component  of  the 
sound  field. 


The  Dynamic  Sea  Surface 

The  representation  (and  justification  therefor  )  which  we  shall  employ 
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hat  bttn  discussed  by  Marsh*  et  al.  Briefly,  we  have  the  isotropic  turbu¬ 
lence  spectrum  H  (K,  <i>).  The  auto  covariance  function  of  surface  elevation 
it  p: 

p(r,/)  »  <Z(x,y,t)Z(x  +  6,  y  +  rj,  t+7^)>  (2) 


where  <  >  represents  an  ensemble  average  and  r*  *  {*  +  r*.  The  turbulence 
spectrum  it  defined  to  be 


H  (K,  (i) ) 


<vo  C>0 


p  (r,  7")  cos  (w  7)  J0  (Kr  )rdrd  7^(3  ) 


This  it  the  ordinary  three  dimensional  power  spectrum;  the  Bessel  function 
So  enters  because  of  assumed  isotropy.  There  are  also  the  marginal  spectra 


oo 


to 


A1  (K) 


H  (K,  u)du>  and  B*  (<*>)=  /  H  (K,u)dK.  (4) 


If  K  and  w  are  functionally  related,  we  have  A2  (K  )dK  n  B1  (w)dw,  and  then 
H  (K,  w)has  values  only  on  the  diagonal  K  x  K  (w  ).  In  this  case  we  speak 
of  dispersion,  and  there  is  a  definite  speed  v  ■  w/K  (w)  for  each  wave 
number.  The  classical  dispersion  relation 


J  •  gK  ♦  y  K*  (5) 

(g  *  acceleration  of  gravity  and  y  the  specific  surface  tension,  about  74cgs 
for  clear  water  )defines  a  minimum  phase  velocity  v  «  (4y  g)  ^4,  which 
corresponds  to  a  wave  frequency  near  14  Hs  and  wave  length  of  about  2  cm. 
Waves  longer  than  this  are  gravity  controlled  and  dispersive: 


A*  (K)  «  p  expt^^/SM/IK* 

Bs  (w)  *  p  g*  exp  (-  2  v*/S*  )/m* 


for  a  wind  speed  S. 


The  dimensionless  constant  p  is  7.  4  x  10*1. 


(6) 
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Short  wi 


o  turbulent*  and 


H  (K*  w)  *  «xp  (•  v1  /2  vi  Y  w/2/K4Vpo  (7) 

where  v  ■  «/K,  vM  ■  v0  (K0/K  f  N,  and  the  numeric  N  is 

-y  2  f/*1  2/*  r  (  n/6)]'^ 

The  treneition  region  near  the  minimum  phase  velocity*  ie  unfortu¬ 
nately  not  well  defined*  although  there  ie  evidence  that  turbulence  extends 
somewhat  into  the  gravity  range. 


Ambient  Pressure 

Marsh1  has  shown  that  the  power  spectrum  of  surface  generated 
noise  is  proportional  to  the  sea  surface  velocity  spectrum  which  is  simply 
JB*  (  w).  In  like  manner*  the  cross  spectra  of  pressures  at  separated 
points  may  be  obtained  from  a  weighted  integral  of  the  turbulence  spectrum* 
with  weight  dependent  upon  the  geometry.  We  are  not  yet  able  to  report 
results  on  these  cross  spectra*  but  will  use  measured  or  postulated  co- 
variance  functions  which  are  consistent  with  the  desired  power  spectra. 

0 

Such  covariance  functions  are  reported  in  the  literature  by  Cron4  *  and 
others.  Modeling  wave  forms  with  prescribed  cross  spectra  will  be  taken 
up  below. 


Propagation 


Thi  dlrtct  arrival  hat  been  de scribed  by  Marsh  et  al5.  The  con¬ 
trolling  equation  is 

iq  +  B/BY)y>  ■  (q  +  8/9  f  Myy . )  +  (qCiy^_  +  C*y_  ...  )  ®XP  6  (B) 

£  7  i  i  i  i  i 

y  is  the  scaled  pressure: 

y  «  2BpR/po 
$  the  scaled  range: 

t  *  ln(R/R0) 

the  local  time: 

T  m  t  m  R/c* 

q,  Ci  and  Ci  are  relaxation  and  viscous  coefficients  characteristic  of  sea 
water;  B  is  determined  by  the  lowest  order  non-linear  terms  in  the  equa¬ 
tion  of  state*  and  po  is  a  reference  pressure.  Integration  of  the  equation 
from  a  prescribed  wave  form  f  (t  )at  reference  range  R0  has  been  described. 
The  same  equation  may  be  used  to  determine  the  coherent  scattered  wave 
if  R  is  taken  to  be  the  total  path  length,  if  the  sign  of  y  is  reversed,  *nd  the 
amplitude  reduced  by  a  factor  1  -  a.  Measurements  and  theory  for  this 
"surface  loss"  have  been  reported?  The  flux  which  is  lost  from  the  coherent 
wave  reappears  in  the  incoherent  wave,  distributed  in  angle  and  frequency. 
The  incoherent  power  spectrum  at  a  point  may  be  obtained  from  a  weighted 
integral  of  the  surface  turbulence  spectrum,  with  weight  dependent  upon 
geometry,  and  in  a  similar  manner,  the  cross  spectra  at  separated  points 
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may  be  obtained.  Here  again*  we  have  no  results  to  quote  on  cross  spectra. 
Two  simple  cases  have  been  investigated*  as  follows.  For  monostatic  re¬ 
verberation*  the  power  spectrum  of  the  reverberation  is 


sR<«> 

M  (t) 


.  OO 

M(t)  I  wj  |e(c^)|  H  (K,  u-ubldufe 


O 


128  Z q  /  c "  cos  8  t7 


K  >  2^  cos  9/c  ;  sin  6  =  2  Z0/  ct 


(9) 


Actually*  since  the  factor  M  depends  upon  time*  this  is  a  pseudo -spectrum. 
Its  nature  and  some  of  its  limitations  are  discussed  by  Faure7.  The  quantity 
E(«)  2 is  the  "energy"  spectrum  of  the  pressure  form  f  (t ) : 


For  short  signals  with  energy  scattered  mainly  by  gravity  waves,  the  pseudo¬ 
power -spectrum  of  the  monostatic  reverberation  is  proportional  to  the  signal 
"energy"  spectrum. 

The  other  example  concerns  forward  scatter.  For  frequencies  be¬ 
low  a  few  kHs,  depending  upon  geometry  and  sea  conditions,  the  incoherent 
component  of  the  field  will  consist  of  upper  and  lower  side  bands  of  spectrum 
proportional  to  the  marginal  surface  frequency  spectrum?  For  an  isotropic 
sea  the  side  bands  will  be  equal  in  intensity*  but  not  if  the  sea  is  directional. 
The  percent  modulation  is*  as  indicated  previously,  determined  by  the  "sur¬ 
face  loss"  in  such  a  way  that  the  total  flux  of  energy  incident  upon  the  surface 
is  recovered  in  the  scattered  field. 
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Simulation 


In  the  sequel*  we  shall  discuss  simulation  without  regard  to  the  par¬ 
ticular  stochastic  elements  involved.  This  will  permit  us  to  update  know¬ 
ledge  of  sea  noise*  propagation*  etc.  and  to  introduce  "real"  data  at  will. 
Although  the  simulation  will  be  electrical*  we  shall  maintain  sound  pres¬ 
sure  units  throughout.  We  deal  with  an  array  of  pressures  (t )  and 
sampled  equivalents 


A 

nm 


p  (m  T ) 
n 


(ID 


which  are  quantised  and  obtained  with  the  sampling  interval  T.  In  our  simu¬ 
lation  to  date  we  have  used  12  bit  quantization  and  intervals  as  small  as  one 
microsecond.  By  way  of  notation*  the  hat  (circumflex)  will  show  that  the 
quantities  is  a  sample  value  of  a  continuous  variable.  On  the  other  hand*  a 
number  derived  (e.  g.  *  by  filtering  )  from  Pnm  will  not  carry  th^  liat.  If  only 
a  single  pressure  function  is  of  interest*  the  first  subscript  will  not  be  written. 


Filtering*  or  spectral  modification  of  p  (t)may  be  simulated  by  a  trans¬ 
versal  digital  filter9.  The  relation 


2  A 

n 


A 

P 

m-n 


+  2  B 

n 


y 


m-n 


(12) 


corresponds  to  filtering  with  a  rational  function  frequency  response.  The 

coefficients  A  *  B  are  determined  by  the  filter  poles  and  zeros.  A  simple 

n  n 

low  pass  filter  requires  2  coefficients;  a  high  pass*  3  and  a  bandpass  4. 
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Phase  Shifter 


The  frequency  response  (X  +  la  )/ (X  -  i(t> )  is  all  pass*  with  phase 
ranging  from  zero  at  low  frequency  to  r  at  high.  A  corresponding  digital 
filter  relation  is 


(1  +  x)y 


m 


( 1  -  X )  +  X  p 


m 


m- 1 


xym 

m  “ 


(  13) 


By  cascading*  and  series  parallel  arrangements  of  several  basic 
phase  shifters*  any  desired  phase  shift  frequency  characteristic  can  be 
approximated*  and  then  multi  channel  wave  forms  with  prescribed  cross 
spectra  can  be  simulated. 


Random  Number  Generator 

Simulation  of  a  stochastic  variable  starts  with  a  random  number  se¬ 
quence.  The  scheme  we  have  used  is  described  in  the  literature10'  ll.  Let 

R  be  a  random  number*  and  (R  )  the  operation  of  shifting  the  number  in 
n  no 

binary  form  a  given  number  of  (binary )  digits  to  the  left.  Then  if  K  is  a 
constant*  and 

R  *  (KR  )  ,  (14) 

n  n-  i  o 

some  K  numbers  may  be  generated  before  repetition.  The  first  number 
should  be  odd;  otherwise  zero  might  develop*  with  early  termination  of 
the  sequence.  It  is  stated  that  K  should  be  a  large  odd  power  of  five.  Our 
Recomp  III  uses  a  40  bit  word.  K  ■  5sf  is  a  satisfactory  choice*  with  left 
shift  of  12  bits  and  retention  of  28  bits.  Other  choices  of  K  and  shift  have 
been  inferior  in  one  respect  or  another.  Starting  with  R0  s  2 15  f  1,  and 
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scaling/Rj*  1,  Rg  »  •  0.  4700.  Table  1  then  gives  the  distribution  of 
amplitudes  and  aero  crossing  intervals  of  the  first  10,  000  numbers.  It 
may  be  seen  that  the  sequence  is  nearly  uniformly  distributed  between  -  1 
and  +  1,  and  the  intervals  between  zero  crossings  approximate  the  desired 
distribution  P^  *  2  n  *.  A  chi-square  test  on  the  tables  gives  a  confidence 
level  of  95%  for  the  desired  distributions. 


Simulation  of  Sea  Noise 


A  simple  choice  is  "pink11  noise,  F  for  which  we  use  a  low  pass 

m 

filter: 

P  =  cR  +  xP  (15) 

m  m  m  - 1 

Assuming  the  R  to  have  the  auto  covariance 

m 

<  R  R  >  s  or2  6  =  tr2,  m-n 

m  n  mn 

s  o,  m  ^  n, 


the  P  have  the  auto  covariance 
m 


<  P  P  >  s 
m  n 


c2  x 


1  -  xr 


I  m-n* 


(16) 


For  the  noise  reported  here,  v2  *  1/  3.  We  have  used  x  =  0,  995, 

c  *  0.  005,  which  corresponds  to  a  power  spectrum  down  3dB  from  its 

low  frequency  maximum  at  about  800  Hs,  for  a  sampling  interval  of 

lp  sec.  Table  II  gives  the  distribution  of  amplitudes  and  sero  crossing 

intervals  of  the  first  10,  000  P  The  maximum  value  of  P  in  this 

m  m 

interval  was  0.  095,  there  being  7  negative  and  5  positive  values  close 
to  this  number. 
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A  more  elaborate,  and  a  realistic  representation  of  sea  noise  re¬ 
quires  the  simulation  of  surface  waves,  with  amplitude  spectrum  B2  («). 
To  date,  we  have  preapred  a  simulation  covering  the  gravity  region,  and 
one  in  the  short  capillary  region.  For  the  first,  we  use  Equation  (6). 


for  large  N.  The  term  can  be  approximated  by  low  and  high  pass 
filters.  Figure  1  gives  an  example,  with  formulae,  using  N  =  40,  for  a 
total  of  48  coefficients  in  the  equivalent  transversal  filter.  In  the  capillary 
region.  Equation  (7  )  leads  to 

B 2  (w)  -  w  7s  (17) 

This  can  be  approximated  by  six  coefficients  as  shown  in  Figure  2. 


Simulation  of  Reverberation 

If  the  approximations  leading  to  Equation  (9  )and  discussed  following 
Equation  (10) are  accepted,  the  reverberation  spectrum  (pseudo) is  propor¬ 
tional  to  the  signal  energy  spectrum.  ”or  an  explosive  source,  the  energy 
spectrum  of  the  main  part  of  the  wave  of  detonation  is  proportional  to 
l/  (4r2f02+c/).  The  cutoff  frequency  f0  is  determined  by  the  explosive 
yield,  hydrostatic  pressure  and,  to  a  small  extent,  range,  and  typically  is 
the  order  of  one  kHz  for  yields  of  a  few  pounds.  Thus,  the  pink  noise  clis  - 
cussed  above  is  also  a  simulation  for  reverberation.  Figure  3  displays  the 
wave  from  short  portions  of  a  reverberation  record  of  3.  8#  of  TNT,  sampled 
at  a  rate  of  40  kHz.  Figure  4  gives  the  distribution  of  amplitude?  of  2nr>0 
sample  points,  and  also  the  distribution  of  pink  noise  with  f0  =  80.1  Hz  and 
the  same  sampling  rate.  The  amplitude  distributions  are  similar,  and  the 
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zero  crossings  are  fairly  alike,  as  may  better  be  seen  in  Figure  5,  which 
shows  the  probability  distributions  of  the  samples.  The  amplitude  distribu¬ 
tion  of  the  noise  is  easily  tailored  by  adding  together  a  few  independent 
random  numbers.  The  law  of  large  numbers  works  rapidly,  and  10  or  so 
numbers  added  together  provide  an  excellent  normal  distribution. 


Summary 

We  have  sketched  here  a  method  of  similating  underwater  sound,  em¬ 
phasizing  the  influence  of  the  dynamic  sea  surface  and  the  use  of  digital 
wave  form  synthesis.  Although  the  methodology  and  results  presented  rep¬ 
resent  hardly  more  than  a  beginning,  it  is  believed  that  known  properties 
of  real  sounds  can  be  accurately  simulated.  Modern  computers  are  so  fast 
that  real  time  simulation  is  now  a  reality,  and  it  is  inevitable  that  digital 
filtering  will  be  of  Increasing  value  for  on-line  processing. 
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Table  I.  Distributions  of  10,000  Random  Numbers  (Ks5,T ) 

Amplitudes 


-  1.  0 

To 

-  0.  9 

470 

-0.9 

To 

-  0.  8 

517 

-  0.  8 

To 

-  0.  7 

498 

.  0.  7 

To 

-  0.  6 

508 

-  0.  6 

To 

.  0.  5 

517 

-  0.  5 

To 

-  0.  4 

481 

-0.4 

To 

-  0.  3 

515 

.0.3 

To 

.  0.  2 

506 

-  0.  2 

To 

-  0.  1 

511 

-  0.  1 

To 

-  0.  0 

494 

.  0.  0 

To 

0.  1 

493 

0.  1 

To 

0.  2 

493 

0.  2 

To 

0.  3 

540 

0.  3 

To 

0.  4 

489 

0.4 

To 

0.  5 

513 

0.  5 

To 

0.  6 

498 

0.  6 

To 

0.  7 

506 

0.  7 

To 

0.8 

494 

0.8 

To 

0.9 

467 

0.  9 

To 

1.  0 

490 

Zero  Crossing  Intervals 

Interval 

Positive 

Negative 

1 

1273 

1235 

2 

625 

690 

3 

328 

300 

4 

150 

167  ; 

5 

97 

68  ' 

6 

31 

33 

7 

15 

17 

8 

4 

10 

9 

4 

6 

10 

2 

2 

11 

0 

1 

12 

1 

1 

13  and  more 

0 

0 

i 

-t 
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Table  II  Distributions  of  10,  000  values  of  Pink  Noise  (K»515 ) 


Amplitude  Number 


0.  10 

To 

0.  09 

7 

0.  09 

To 

0.  08 

39 

0.  08 

To 

0.  07 

189 

0.  07 

To 

0.  06 

248 

0.  06 

To 

0.  05 

477 

0.  05 

To 

0.  04 

756 

0.  04 

To 

0.  03 

819 

0.  03 

To 

0.  02 

880 

0.  02 

To 

0.  01 

998 

0.  01 

To 

0.  00 

1324 

0.  00 

To 

0.  01 

1263 

0.  01 

To 

0.  02 

981 

0.  02 

To 

0.  03 

598 

0.  03 

To 

0.  04 

505 

0.  04 

To 

0.  05 

289 

0.  05 

To 

0.  06 

326 

0.  06 

To 

0.  07 

193 

0.  07 

To 

0.  08 

83 

0.  08 

To 

0.  09 

20 

0.  09 

To 

0.  10 
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Zero  Crossing  Interval 

No. 

Positive 

Negative 

1 

54 

65 

2 

22 

17 

3 

13 

13 

4 

12 

9 

5 

6 

4 

6 

5 

5 

7 

2 

2 

8 

3 

2 

9 

4 

2 

10 

1 

4 

11 

1 

2 

12 

0 

0 

13 

1 

0 

14 

2 

3 

15 

2 

1 

16 

1 

0 

17 

5 

4 

18 

1 

1 

19 

3 

0 

20 

0 

3 

13 


Zero  Crossing  Interval 


No. 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

>  30 


Positive 

3 

1 

1 

1 

1 

1 

0 

0 

0 

0 

47 


largest  interval  1138 


Negative 

0 

0 

0 

1 

2 

1 

0 

0 

0 
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APPROXIMATION  TO  GRAVITY  WAVE  SPECTRUM. 

LINE  IS  DESIRED  FORM;  POINTS  ARE  FROM  APPROXIMATING  FORMULA. 
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TIME  (msec) 

WAVE  FORM  OF  12.5  ms  OF  SURFACE  REVERBERATION  FROM  3.8  POUNDS  OF  TNT. 
SAMPLES  WERE  TAKEN  SHORTLY  AFTER  THE  BEGINNING  OF  THE  SURFACE  RETURN 
(TOP  TRACE)  AND  SOMEWHAT  LATER,  WELL  DOWN  ON  THE  DECAY  CURVE  (BOTTOM  TRACE). 
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AMPLITUDE 

DISTRIBUTIONS  OF  AMPLITUDES  OF  REVERBERATION  AND  PINK  NOISE. 
AMPLITUDE  SCALE  FOR  "LATE"  REVERBERATION  IS  TEN  TIMES  THAT  OF  EARLY 
REVERBERATION  AND  PINK  NOISE. 
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Figures 


1.  Approximation  to  Gravity  Wave  Spectrum. 

Line  ie  desired  form;  points  are  from  approximating  formula. 

2.  Approximation  to  Capillary  Wave  Spectrum. 

Line  is  desired  form;  points  are  from  approximating  formula. 

3.  Wave  form  of  12.  5  ms  of  surface  reverberation  from  3.  8 
pounds  of  TNT.  Samples  were  taken  shortly  after  the  be¬ 
ginning  of  the  surface  return  (top  trace)  and  somewhat  later, 
well  down  on  the  decay  curve  (bottom  trace  ) . 

4.  Distributions  of  amplitudes  of  reverberation  and  pink  noise. 
Amplitude  scale  for  Mlatd' reverberation  is  ten  times  that  of 
early  reverberation  and  pink  noise. 

5.  Distribution  of  sero  crossings  of  reverberation  and  pink  noise. 
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Spectra  of  The  Dynamic  Sea  Surface 


H.  W.  Marsh*  and  R.  H.  Mellen* 


Abstract 


A  unified  model  of  the  equilibrium  surface,  including 
turbulence,  is  proposed.  The  me  cl  is  based  upon  inter¬ 
pretations  of  a  variety  of  measurements,  including  physi¬ 
cal  (wave  staff);  optical  glitter  and  slope;  radar  and 
sonar  scatter  and  acoustic  ambient  noise. 
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SPECTRA  OF  THE  DYNAMIC  SEA  SURFACE 


Introduction 


Experimental  measurements  of  the  statistical  behavior  of  the  wind- 
driven  sea  surface  have  proven  surprisingly  difficult  to  reconcile  to  a  uni¬ 
fied  model  A  list  of  some  of  the  various  methods  and  sources  is  given  in 
Table  L 


TABLE  I 


Category 


Type  of  Measurement 


Source 


Physical 


Optical 


Optical 


Wave  height  vs  frequency  by  wave 
staff  measurements,  etc. 


Slope/angle  distribution  by  sur¬ 
face  glitter  photography 


Slope  spectrum  vs  frequency  mea¬ 
surements  by  refraction  of  light 


Burling  (  19  55) 
Neumann  (  19  53) 
Pierson  (  1955) 


Cox  (  19  56) 


Cox  (  19  58) 


Radio 


Radar  scatter  strength  vs  wave  — 
length  and  angle,  Doppler  Fre¬ 
quency  spectra 


Kerr  ( 19  5 1) 

MacDonald  (1962) 
Wiltsie  (  19  57) 


Underwater  Sound  Sonar  scatter  strength,  Doppler 

spectra,  ambient  noise  spectra 


Marsh  (1963a,  b) 

Marsh  (1964) 

Mellen  (1964) 


Sound  in  Air 


Scatter  strength,  Doppler  spectra  Liebermann  ( 1963  ) 
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The  main  difficulty  is  that  each  case  requires  appropriate  theoreti¬ 
cal  interpretation.  Thus,  electromagnetic  and  acoustic  scattering  at  "fi¬ 
nite"  wavelengths,  optical  refraction  and  reflection,  ambient  noise,  and 
Doppler  theories  are  all  involved  in  relating  these  data  to  the  dynamics 
of  the  surface.  We  now  believe  that  Marsh*s  ( 19 6 1 )  first  order  scattering 
theory,  the  surface  noise  radiation  theory  (Marsh  1963a)  and  Kuo's  (1965) 
hydrodynamic  theory  of  the  turbulent  surface,  together  with  the  others, 
give  a  reasonably  complete  picture. 

According  to  Kuo,  both  turbulent  surface  flow  and  wave  motion  are 
expected  to  grow  under  the  influence  of  the  turbulent  wind.  The  equilibrium 
excitation  levels  are  known  only  by  experiment  since  the  theory  does  not 
consider  the  saturated  state.  We  find  by  sonar  and  radar  scattering  mea¬ 
surements  that  the  height  spectrum  is  asymptotic  to  K"1,  independent  of 
wind  speed  (figure  l)  .  This  gives  a  scattering  strength  that  is  also  inde¬ 
pendent  of  wavelength  (figure  2)  . 

♦ 

From  Doppler  measurements  at  short  wavelengths  (figure  3)jwe 
also  find  that  the  associated  frequency  or  horizontal  velocity  spectra  are 
not  related  to  the  wavenumber  spectrum  by  the  gravity-surface  tension 
dispersion  formula.  Instead,  the  surface  appears  to  act  more  like  a  tur¬ 
bulent  boundary  layer.  At  very  long  wavelengths,  however,  radio  mea¬ 
surements  do  show  Doppler  frequencies  associated  with  the  resonant  gravi¬ 
ty  wavelengths.  Also,  from  forward  scatter  Doppler  calculations  (figure 
4)  ,  we  find  a  similar  situation  since  the  Doppler  spectrum  is  dominated 
by  the  longest  wavelengths,  and,  indeed  is  the  weighted  surface  frequency 
spectrum. 


♦ 

See  figures  and  discussion  and  conclusion  below. 
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In  the  gravity  region  according  to  surface  wave  theory,  the  K 
wavenumber  spectrum  transforms  to  the  w"1  frequency  spectrum,  which 
agrees  with  results  obtained  directly  from  wave  height  measurements. 

In  the  capillary  region,  however,  the  experimental  results  do  not  agree 
with  wave  theory.  The  frequency  spectrum  is  derived  from  both  Doppler 

- 11/  -7/ 

and  ambient  noise  measurements  and  gives  o>  /3  instead  of  □  /3. 

The  transition  between  gravity  waves  and  turbulent  flow  is  not  known 
with  any  exactness.  In  fact,  the  surface  may  be  truly  wavelike  only  for 
the  longest  existing  wavelengths.  Complete  turbulence  may  be  estimated 
to  begin  somewhere  near  the  group  velocity  minimum  for  which  K  =202m  1 
However,  the  transition  appears  to  involve  most  of  the  gravity  spectrum. 

It  does  not  appear  likely  that  specific  surface  details  such  as 
breaking  waves,  spume,  etc.  ,  will  alter  the  picture  because  all  these  ap¬ 
pear  to  be  included  in  the  statistical  description  of  the  random  dynamic 
surface. 


1 


Theory 


Observations  and  analyses  drawn  from  a  number  of  sources  have 
been  interpreted  to  provide  an  estimate  of  the  wave  number  and  frequency 
spectra  of  the  wind  driven  sea  elevation.  A  brief  summary  is  given  here, 
covering  frequencies  up  to  5  x  104  Hz  and  wave  numbers  to  l04m  1 .  The 
results  are  given  for  a  fully  developed  sea,  in  equilibrium  with  the  local 
wind,  in  terms  of  an  equivalent  isotropic  spectrum. 


The  xy  plane  corresponds  to  the  mean  surface.  The  instantaneous  el¬ 
evation  z  =  £(x,  y,  t)  has  the  auto  covariance  p  =  <  £  (x,  y,  t )  £(x+£,  y+  rj,  t  +  r )  >. 

A  stationary,  homogeneous,  isotropic  statistical  surface  is  postulated,  and 

l  / 

then  p  depends  only  upon  r  and  r  =  (£2  +  rj2  )  .  A  convenient  two  dimen¬ 

sional  spectrum  is 


H  (K,  u>)  ss 


(r,  t  )  cos  (i)  T  J0  (Kr  )  rdr  dT 


which  has  the  inverse  relationship 


.00*00 


P  (r,r  )  a  f  f  H{K,  w)  cos  wr  J0  (Kr)  d  K  d  t* 

J0  •o 


Thus,  the  variance  (mean  square  wave  height) ,  h2  =  p  (o,  o) ,  is  distributed 
over  the  spectrum  H  (K,  w)  ;  there  are  the  marginal  spectra 

AJ(K)  =/*  H(K,u)dw=  -=r-J  p  (r,  o)  J„  (Kr)  rdr 

/a  c.  n  4 


and 


B2  ( w)  =  f  H  ( K,  w)  dK  =  -1-  T®  p  ( o,  t  )  cos  w  r  dT  . 

Jo  n  j  a 
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The  variance  is  given  in  terms  of  either  of  the  marginal  spectra: 

oo 

h2  =/  A2 ( K ) dK  -t  B2  (u)  du  . 

There  is  the  dispersion  equation  for  surface  waves  in  deep  water 
(neglecting  viscosity )  which  relates  w  and  K: 

=  gK  +  yK3 

Here,  g  is  the  acceleration  of  gravity  and  y  the  specific  surface  tension 
(about  74  cgs  for  clean  water) .  When  the  dispersion  relationship  applies, 
H  (K,  u>)  is  relevant  only  on  the  diagonal  corresponding  to  K  =  K0  (o>) 

H(K,u)  =  A2(K)6(K-K0)=  B2  (u)  6(K.K0). 

K0  satisfies  the  dispersion  equation  and  6  is  the  Dirac  function.  In  any 
event,  the  relation 


dw  A2(K) 
dK  =  B2  (u) 


defines  a  functional  relationship  between  K  and  which  we  may  call  the 
mean  dispersion.  If  H  (K,  w)  is  not  diagonal,  the  surface  elevation  is  tur¬ 
bulent  instead  of  wavelike.  In  this  case,  there  is  no  deterministic  dis¬ 
persion  equation,  but  only  mean  dispersion. 


For  dispersive  waves,  the  phase  velocity  is  v  =  w/K  =(g/K  +  yK) 

«  /  i  / 


v, 


1/ 

The  minimum  velocity  v0  =  (4  y  g)  ,  corresponding  toK0  =  (g/y)/2  and 
«b  =  (4g3/y)^4.  Waves  for  which  K  <  K0  are  the  gravity  waves;  those 
for  which  K  >  K0  are  the  capillary  waves. 


V; 
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The  Spectra 


I.  Gravity  region:  In  the  gravity  region,  dispersive  waves  pre¬ 
dominate.  The  spectra  are  of  the  form 

A2  (K)  =  p  exp  (-  2  v2 /s2  )  /2K3 
B2  ( <0  )  =  p  g2  exp  (-2v2/s2)/u5 

for  a  wind  speed  s.  The  dimensionless  constant  (3  is  7.  4x10  3,  ac¬ 
cording  to  Burling  (  1955  ).  These  spectra  are  in  general  agreement  with 
the  results  of  Pierson  (  1964). 


II.  Capillary  region:  In  the  capillary  region,  turbulence  pre¬ 
dominates.  The  marginal  spectra  are  found  to  be 

A2(K)  *  p /2KJ 

B2  (  u  )  =  Pg2/«o^  u  ^ 


Note  that  the  Burling  spectrum  persists;  the  frequency  spectrum  has 
changed  form  in  accordance  with  the  shift  from  gravity  to  surface  tension 
control.  The  entire  spectrum  H  ( K,  w)  is  displayed  in  cross  section  in 
Figure  2,  (see  discussion  and  conclusions  below).  These  may  be  ap¬ 
proximated  by  the  formula 


H(K,  w)  =  exp  (  -  v2  /2  v  200  )  (3/  ^/tt/2  K4  v0o 


1/ 

where  v  =  w/K  and  Voo  =  Vo  ( K0  /K  )  N.  The  numeric  N  is  equal  to 

5/6_  7s 


/16  (2)/6  T  (11/6)]  8  . 
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Discussion  and  Conclusions 


Figure  1  is  an  idealized  height-wavelength  spectrum  for  two  wind 
speeds,  5  and  10  ms  ,  illustrating  the  asymptotic  K  dependence  in 
the  saturated  region  and  wind  speed  dependence  of  cutoff. 

Sonar  and  radar  measurements  of  the  sea  surface  fcackscatter  ap¬ 
pear  to  show'  independence  of  wavelength  and  the  predicted  angular  de- 
pendence  of  scattering  strength  consistent  with  the  asymptotic  K  height 
spectrum.  Figure  2  illustrates  the  wavelength  independence  as  well  as 
the  angular  dependence  of  scattering  strength  for  a  rather  wide  range  of 
wavelengths 

The  Doppler  spectra  for  several  wavelengths  reported  by  various 
investigators  are  summarised  in  Table  II.  In  Figure  3}the  data  are  re¬ 
ferred  to  the  equivalent  surface  wavenumber  K  and  the  frequencies  Aw 
are  referred  to  the  apparent  horizontal  surface  velocity  by  V  =  Aw/K. 

The  normalized  Doppler  spectra  are  shown  in  dB  relative  to  the  maximum 
which  is  10  dB  above  the  baseline.  Where  two  spectra  are  shown  sym¬ 
metrically  about  zero  velocity  (B,  C,  E,  F,  I,  J,  K) ,  the  original  data  did  not 
distinguish  the  sidebands.  The  normalizing  procedure  is  based  on  the  work 
of  Liebermann  (  19  63)  and  Marsh  (  19  63  b) ,  using  first  order  resonance  theory, 
The  dotted  curves  are  ^  V,  the  classical  dispersion  relation.  Actual  spec¬ 
tral  distributions  would  be  governed  by  the  radiation  beam  pattern.  The 
lack  of  agreement  between  the  experimental  data  (A,  C,  D,  E,  F,  G,  H)  and 
classical  theory  is  believed  due  to  surface  turbulence. 

In  Figure  4  is  shown  the  predicted  forward  scatter  Doppler  spectrum 
compared  with  the  actual  surface  frequency  spectrum.  The  similarity  is 
obvious  and  consistent  with  other  forward  scatter  Doppler  measurements. 
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List  of  Illustrations 


Figure  No. 


1.  Sea  surface  elevation  wavenumber  spectrum  for 

wind  speeds  of  5  and  10  meters  per  second. 


2.  Backscattering  coefficients  per  unit  area 


3.  Surface  wave  velocity  from  backscatter  data 


4.  Doppler  spectrum  of  forward  scatter 
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WAVENUMBER  ,  K  (rtf') 

SEA  SURFACE  ELEVATION  WAVENUMBER  SPECTRUM  FOR  WIND  SPEEDS  OF  5  AND  10  METERS/SEC. 

FIGURE  I 


BACK  SCATTERING  COEFFICIENTS  PER  UNIT  AREA 


EQUIVALENT  SURFACE  WAVE  NUMBER.  radions/m 
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APPARENT  SURFACE  VELOCITY,  ( m  s'1 ) 


SURFACE  WVE  VELOCITY  FROM  BACKSCATTER  DATA 


FIGURE  3 
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DOPPLER  FREQUENCY,  Hz 

DOPPLER  SPECTRUM  OF  FORWARD  SCATTER 


System  Parameters  for  Ultrasonic  Explosive  Echo-Ranging 


H.  W.  Marsh  * 


Abstract 


Calculations  are  presented  illustrating  the  use  of  high 
frequencies  (up  to  40  kHz)  for  ranging  on  postulated  tar¬ 
gets.  The  high  frequencies  propagate  with  acceptable  los 
because  of  the  finite  amplitude  of  the  explosive  source. 
Comparatively  high  directional  discrimination  against 
reverberation  and  noise  can  thus  be  achieved  with  a  re¬ 
ceiver  of  small  dimensions. 
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SYSTEM  PARAMETERS  FOR  ULTRASONIC  EXPLOSIVE  ECHO-RANGING 


INTRODUCTION 


The  system  concept  entails  a  3.8  lb.  Sofar  charge,  (other  sizes 
can  be  used),  as  a  source  and  a  directional  receiver  operating  at  ultra¬ 
sonic  frequencies,  (  several  tens  of  kilocylces  per  second  ).  The  novel 
feature  of  the  system  is  its  high  frequency,  compared  with  current  prac¬ 
tice.  Because  of  the  high  frequency,  directivity  is  obtainable  with  a 
small  receiver,  providing  good  resolution  for  discrimination  against  re¬ 
verberation,  and  ambient  noise.  Tests1-2,  at  sea  have  proven  that  the 
desired  frequencies  propagate  with  much  lower  loss  than  is  true  of  small 
amplitude  sounds,  and  indicate  that  very  high  signal  to  interference  ratios 
are  achievable  for  targets  at  ranges  of  interest. 


No  attempt  has  yet  been  made  to  optimize  system  parameters. 
Also,  full  scale  tests  with  an  actual  target  using  an  adequate  receiving 
system  have  not  been  made.  However,  there  is  sufficient  factual  infor¬ 
mation  now  available  to  permit  estimation  of  the  performance  of  postula- 
ted  systems. 


Estimates  of  performance  for  the  system  shown  in  block  diagram, 

(Figure  1),  has  been  made  using  the  simulation  procedures  described  by 
Marsh3 . 
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The  performance  criterion  used  is  peak  signal  response  (A  in  the 
Figure),  and  peak  interference  (B  in  the  Figure  )  for  a  specified  obser¬ 
vation  time. 


SIGNAL 


The  signal  input  to  the  receiver  is  a  recording  of  the  direct  arrival 
from  a  3.8#  Sofar  bomb  detonated  at  a  range  of  6150  meters  (the  leading 
edge  is  shown  as  I,  Figure  9  of  Reference  1  ).  The  shape  of  the  main 
arrival  from  an  actual  target  at  a  range  of  4000  meters  is  believed  to  be 
represented  by  this  form.  The  rise  time  is  kept  small  by  finite  ampli¬ 
tude  effects  during  the  4000  meter  path  to  the  target,  and  then  increased 
normally  during  the  return  trip.  The  rise  times  tj  out  and  t2  back  will 
tend  to  add  in  quadrature  (  see  description  following  Equation  52  of  Ref¬ 
erence  1)  to  produce  a  net  rise  time  t3  at  the  receiver  of  (ti2  +  t22 )  =  t3 . 

Peak  pressure  (broad  band  )at  the  target  is  100  db//|ibar  (Figure  10  of 
Reference  1).  Loss  in  peak  amplitude  on  return  trip  is  estimated  to  be 
,77  dB  (Figures  6  and  10  of  Reference  2).  The  effective  target  strength 
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is  more  difficult  to  estimate.  Figures  2-4  show  broad  band  target  echos 
from  a  beam  aspect  submarine,  at  1500  meters  from  l/2#  TNT  source. 

The  principle  echo,  which  is  the  main  event  in  Figure  2  and  can  be  identi¬ 
fied  by  arrows  in  Figures  3  and  4,  has  an  effective  target  strength  of  25dR. 
The  shape  of  the  echo  is  a  close  approximation  to  that  of  the  first  positive 
phase  of  the  direct  arrival  from  the  source,  which  suggests  that  the  target 
reflection  can  be  taken  to  be  distortionless,  but  with  reduction  of  amplitude 
for  aspect  dependence.  Model  studies  reported  by  Rubega4  indicate  that 
the  target  strength  for  a  quarter  aspect  fleet  hull  should  be  about  +7  dB 
for  the  short  pulses  contemplated  here.  With  these  assumptions  the  prin¬ 
ciple  signal  arrival  at  the  receiver  has  an  amplitude  of  30dB//jjtbar. 


DIRECTIONAL  RECEIVER 


An  18"  diameter  spherical  lens  filled  with  CC14  produces  directivity 
patterns  shown  in  Figure  5.  These  were  measured  using  a  "Homemade" 
lens  fabricated  from  hand  laminated  fiberglass.  The  patterns  correspond 
closely  to  the  theory  and  measurements  of  Toulis5.  The  gain  shown  is 
relative  to  OdB  for  an  omnidirectional  hydrophone.  This  gain,  at  maximum 
response  can  be  approximated  by  a  high  pass  filter  with  cut-off  frequency 
of  3500  Hz  and  gain  of  25  dB  at  infinite  frequency*. 


Tests  recently  completed  show  that  the  gain  is  essentially  constant  at 
about  +25  dB  from  20  kHz  to  at  least  100  kHz. 
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FILTERS 


The  high  pass  filter  is  taken  to  have  a  cut-off  frequency  of  800  Hz, 
and  the  band  pass  to  have  a  center  frequency  of  40  kHz  with  Q  of  2.  5 
(band  width  =  16  kHz).  The  high  pass  is  perhaps  unnecessary,  but  is 
used  in  the  simulation  process  to  make  "pink  noise",  or  in  other  words, 
to  flatten  the  ambient  noise  spectrum  below  800  Hz.  The  filter  Q  appears 
optimum  for  the  stated  center  frequency,  but  the  latter  is  only  postulated. 
Finally,  the  phase  equalizer  corrects  for  phase  shift  in  the  pass  band, 
producing  some  signal  enhancement  with  no  effect  upon  noise. 


NOISE  SIMULATION 


Following  Reference  3,  an  isotropic  noise  model  is  used,  with  a 
spectrum  of  the  form  ( 1  +  F2  T2 )  \  3  dB  down  at  800  Hz.  This  corre¬ 

sponds  either  to  typical  ambient  noise  or  to  reverberation  with  a  cross 
section  independent  of  frequency  (for  the  source  assumed  here). 


RESULTS  OF  SIMULATION 


Table  I  shows  the  effect  of  each  component  of  the  receiving  system 
upon  signal  and  noise. 
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TABLE  I 


SIGNAL  NOISE  (  Sea  State  2  ) 


Input  level,  peak 

+  30.  0  db/jibar 

m 

Output  from  lens 

+  55.  0  db/jibar 

- 

Output  from  HP  filter 

+  54.  0  db/fibar 

+  1.0  db/fibar  * 

Output  from  BP  filter 

+  41.  0  db/^bar 

-  20.  0  db/pbar 

Output  from  0  equal. 

+  43.  0  db/^ibar 

-  20.  0  db/fibar 

The  input  peak  signal  to  peak  noise  after  HP  filter  is  +28  dB;  the 
output  ra':io  is  +  63  dB.  This  entails  an  overall  processing  gain  of  3  5  dB, 
of  which  '25  uB  car  be  attributed  to  the  directional  receiver  and  10  dB  to 
filtering. 

These  results  can  be  interpreted  in  terms  of  reverberation  as  the 
limiting  inte  .  r  •  ■  uce.  A  typical  explosive  reverberation  record  is  shown 
in  Figure  6.  The  shot  was  detonated  at  a  depth  of  4000  feet  and  recorded 
on  a  nearby  omnidirectional  hydrophone.  This  particular  trace  was  ob¬ 
tained  by  filtering  with  a  16-20  kHz  bandpass  filter,  but  the  trace  also 
typifies  results  to  be  expected  for  the  postulated  32-48  kHz  filter  and 
the  pressure  level  scale  has  been  adjusted  to  represent  the  latter  filter. 

The  simulated  target  return  can  be  seen  in  proper  level  and  time  relation¬ 
ship  to  the  trace.  The  limiting  interference  in  this  case  is  bottom  rever¬ 
beration.  (Note,  surface -bottom -surface  arrival  just  before  target  return). 


rms  level  is  -  9.  7  dB//jxbar.  Peak  levels  are  for  1.  25  sec  processing  time. 
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The  water  depth  was  1800  fathoms  -  a  more  favorable  geometry  could  have 
been  achieved  with  a  deeper  source  and  receiver. 


p 
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RESULTS  FOR  OTHER  RANGES 


The  situation  illustrated  in  Figure  6  is  believed  to  be  a  realistic 
example,  for  the  chosen  range  of  4000  meters.  The  signal  to  reverbera¬ 
tion  level  difference  there  is  -I-  25  dB;  the  signal  to  noise  level  difference 
was  shown  to  be  +  63  dB.  As  the  range  opens  the  signal  to  reverberation 
level  difference  will  diminish  it  the  rate  of  3  dB  per  distance  doubled, 
(boundary  reverberation)  until  volume  reverberation  becomes  dominant, 
which  would  be  at  about  50  km  for  the  situation  illustrated.  At  the  3  dB 
rate,  the  signal  would  still  be  15  dB  over  reverberation  at  40  km  (for  the 
+7  dB  target).  Signal  to  noise  level  d.  .erence  will  decay  more  rapidly. 

The  centre  frequency  should  be  decreased  by  about  l/2  for  each  four  fold 
increase  in  range,  to  compensate  for  increasing  absorption.  This  will 
produce  a  loss  of  6  dB  due  to  increased  noise  level  and  a  loss  of  6  dB  due 
to  decreased  direction.  Thus  optimized,  for  a  four  fold  increase  in  range, 
the  signal  to  noise  level  difference  would  decrease  by  36  dB  from  +63  to 
+27  dB;  at  32  km  the  S/N  would  be  +9  dB.  With  the  reduced  frequency,  the 
noise  and  reverberation  levels  would  be  about  the  same,  at  the  32  km  range. 
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FIGURES 


1.  System  Block  Diagram 

2.  Broad  Band  Echo  (50  psec/cm) 

3.  Broad  Band  Echo  (  500  psec/cm) 

4.  Broad  Band  Echo  (400  m  sec/cm) 

5.  Beam  Patterns  of  18n  Spherical  lens 

6.  Typical  Explosive  Reverberation  With  Expected  Target  Echo 
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BROAD  BAND  ECHO  (50^SEC/tM) 

FIGURE  2 
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FIGURE  3 


BEAM  PATTERN  OF  18"  SPHERICAL  LENS 

FIGURE  5 
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